Dynamically Stabilized Bright Solitons in a Two-Dimensional Bose-Einstein 

Condensate 
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We demonstrate that matter-wave bright solitons can be stabilized in 2D free space by causing the 
strength of interactions to oscillate rapidly between repulsive and attractive by using, e.g., Feshbach 
resonance. 
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A pendulum with an oscillating pivot can have a stable 
configuration with its bob situated above the pivot (the 
inverted pendulum This well-known but surprising 
phenomenon is attributed to a net stabilizing force pro- 
duced by alternating stabilizing and destabilizing forces 
at a frequency much faster than the natural oscillation 
frequency of the fixed-pivot pendulum. Similar physics 
has been employed to focus beams of charged particles 
in a synchrotron (alternating-gradient focusing j2|) and 
to trap ions in the Paul trap S. In this Letter, we 
present a novel application of such a stabilizing mech- 
anism to producing bright solitons in a two-dimensional 
(2D) Bose-Einstein condensate (BEC). By a 'bright' soli- 
ton we mean a stable, solitary wave whose density is 
greater than the background one. 

In ID, when the interatomic interaction is attractive, 
bright solitons can be stabilized even without a trapping 
potential {§[§]. In 2D or 3D free space, however, the 
kinetic pressure and the attractive force cannot balance, 
and the condensate always collapses or expands. This 
can be understood by the following simple argument. 
When the characteristic size of the condensate is R, the 
kinetic and interaction energies are proportional to R~ 2 
and —R~ d in d dimensions, and an effective potential for 
R is the sum of these energies. The effective potential, 
therefore, has a minimum only for d — 1. Here we present 
a novel method to stabilize solitons in 2D by causing the 
interaction to oscillate rapidly using, e.g., Feshbach res- 
onance f?J. 

The system considered here is a BEC confined in a 
quasi-2D axisymmetric trap ||, where the axial confine- 
ment energy huj z is much larger than the radial confine- 
ment and interaction energies. We assume that the con- 
densate wave function ip is always in the ground state of 
the harmonic potential with respect to the z direction, 
and that the dynamics are effectively 2D in the x-y plane; 
we will justify these assumptions later. We let the radial 
confinement frequency cv±(t) and the s-wave scattering 
length a(t) vary in time. The system is then described 
by the Gross-Pitaevskii (GP) equation: 
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r 2 iP + g(t)M 2 iP, (1) 



scribes the strength of interactions. In Eq. (Q), length, 
time, frequency, and tp are measured in units of do = 
[h/(mu; ± o)} 1/2 , wJq, wxo, and iV 1 / 2 /^, where u ±0 = 
wl(0). 

We make the strength of interaction oscillate rapidly 
at frequency f2. The oscillation of the scattering length 
with fl ~ uij_ is studied in Ref. in a different context. 
Experimentally, the speed of the change in the strength 
of interaction using Feshbach resonance is limited by that 
of the applied magnetic field ~ 1 G//iS jlO| , which is suf- 
ficient for our purpose. In order to avoid nonadiabatic 
disturbances that destabilize a soliton state, we gradually 
switch on interaction and simultaneously turn off the ra- 
dial confinement potential as follows 



g(t) = f(t)(jgo + giBmSU), 
>i(t) = l-f(t), 



where f(t) is a ramp function 
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t/T {0<t<T) 
1 (t > T). 
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where r 2 = x 2 + y 2 and g(t) = (&irmuj z /K) 1 / 2 Na(t) de- 



The interaction and the potential can be changed in- 
dependently using an optical trap and magnetic-field- 
induced Feshbach resonance. 

We numerically solve the GP equation (|l|) with time- 
dependent parameters (U)-(0) using the Crank-Nicholson 
scheme Jll| . The initial state is assumed to be the nonin- 
teracting ground state in the presence of the radial con- 
finement potential with luj_q. We gradually increase the 
strength of interaction and switch off the trap according 
to the linear ramp (|j) with T — 20. Figure [j] demon- 
strates the dynamic stabilization of a soliton in 2D, where 
the parameters used are go = — 2tt, g\ = 87r, and = 40 
[Fig. |] (a)] and tt = 30 [Fig. (b)]. Since the value 
of go = — 27r exceeds the critical strength of attractive 
interaction for the collapse g cr ~ —5.8 [ fl2| , the conden- 
sate would collapse if the oscillating term g\ sin fit were 
absent. Even after the radial confinement potential van- 
ishes (t > 20), the soliton shape is maintained as shown 
in the insets of Fig. 0. The plots of the peak density 
\ip (r = 0)| 2 and the monopole moment (r) = f r\ip\ 2 dr 
rapidly oscillate at frequency ft. We note here that the 
dynamics of BEC can be separated into two parts: a 
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FIG. 1: Time evolution of the peak density \ip(r = 0)j 2 
and the monopole moment (r) — J r\ip\ 2 dv. The interaction 
g(t) — —2n + 87rsinnt is switched on and the radial confine- 
mentnotential is switched off from t — to 20 according to 
Eq. (§), where Q = 40 in (a) and fi = 30 in (b). The plots 
oscillate rapidly at Q (which is beyond the resolution limit of 
the graph), and the vertical width of the plots represents the 
amplitude of the rapid oscillations. The insets show the den- 
sity profiles |i/>(r)| 2 at t = 50 (solid curves) and t = (dashed 
curves) . 



rapidly oscillating part with small amplitude and a slow, 
smoothly varying part. After the rapid oscillations are 
averaged out, the width and the peak density are fairly 
constant as seen in Fig. [I] (a). However, when the nona- 
diabaticity of ramp function f(t) is not negligible, the 
lowest breathing mode (but no other radial modes) is 
excited as illustrated in Fig. |l| (b), where the breathing 
mode is caused by an effective potential due to the os- 
cillating interaction [see Eqs. (|9|)-(|l"l"|)]. The amplitude 
of the excitation decreases with increasing ramp time T. 
Such excitations can be suppressed by choosing the ap- 
propriate ramp functions for g(t) and u>±(t). We have 
confirmed that no symmetry-breaking mode grows in the 



azimuthal direction. The system is therefore free from 
such modulational instabilities. 

We checked the validity of the quasi-2D approxima- 
tion by numerically integrating the GP equation in 3D 
axisymmctric systems with large asymmetry parameters 
u) z /u)j_ ~ 50, and found that soliton stabilization occurs 
without axial modes being excited. Our prediction of dy- 
namically stabilizing bright solitons in 2D is, therefore, 
experimentally feasible with an oblate trap as realized by 
the MIT group §. 

To understand the behavior in Fig. [I], we employ the 
variational method with a Gaussian wave function g [IJ, 
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where R{t) is the variational parameter characterizing 
the size of the condensate, and the imaginary term in the 
exponent describes mass current. Substituting Eq. (j^) 
into the action that derives the GP equation ([l]) and 
minimizing the action with respect to R(t), we obtain 
the equation of motion for R(t) as 



R(t) 



1 go + gi sin fit 



R3(t) 



2irR 3 (t) 



(6) 



where we set f(t) = 1. We separate R(t) into the slowly 
varying part Ro(t) and the rapidly oscillating part p(t) 
as R(t) = Ro(t) + p{t) §. When fi > 1, p(t) becomes 
of the order of S!~ 2 , and we keep the terms of the order 
of up to fl~ 2 in the following analysis. Substituting this 
R(t) into Eq. (^), we obtain the equations of motion for 
the rapidly and slowly varying parts as 
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Ro(t) 



91 



27TR 3 (t) 



sin Qt, 



9o 



2tt 



35i 



2itRl{t) 2vri?4(t) 



p(t) sin Sit, 
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where the overline indicates the time average of 
the rapid oscillation. Equation ([?]) yields p(t) — 
-gi/[27rn 2 i?j^)]sinf^. Substituting this into Eq. (§), 
we obtain the equation of motion for the slowly varying 
part as 
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dU(Rp) 
dRo 



(9) 
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The minimum of the effective potential U is attained at 
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and the frequency of small oscillations (breathing mode) 
around the minimum is given by 
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(11) 
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FIG. 2: (a) (gri/fi) 1 / 2 dependence of the monopole moment 
(r) and (b) il/gi dependence of the breathing- mode frequency 
uibr, where go — —2-k. The circles are obtained by varying g\ 
for SI = 30, and the squares by varying Q for g± — 8n. The 
dashed lines show 1.06(gi/Q,) 1 ^ 2 in (a) and 0.32f2/gi in (b), 
which are determined to best fit the plots. 



The Gaussian approximation thus indicates that a soliton 
is stable for go < —2tt. 

The physical mechanism of soliton stabilization can be 
understood from the above discussion. In the inverted 
pendulum, the interplay between the micromotion of the 
bob and the force gradient (i.e., stronger oscillating force 
for larger deviation from the equilibrium position) pro- 
duces a pseudopotential. Since the pseudopotential is 
proportional to the square of the amplitude of the os- 
cillating force a potential barrier is formed around 
the equilibrium position, thereby preventing the pendu- 
lum from swinging down. Such a mechanism also sta- 
bilizes a BEC in a double-well potential with oscillating 
interaction [|l5) and an optical beam propagating in non- 
linear medium with an alternating nonlinearity p6|. In 
the present case, the oscillating "force" for R is given 
by gi sin Jlt/(27ri? 3 ), which becomes larger for smaller R. 
This force gradient and the micromotion of R produce a 
pseudopotential proportional to the square of the "force" 
oc i?~ 6 , which prevents the system from collapsing by 
counteracting the —R~ 2 term that describes the attrac- 
tive interaction. 

In order to stabilize the soliton, \g \ must exceed the 
critical value of collapse, |g C r| — 5.8, which is smaller than 
that of the Gaussian approximation, |g cr | = 2n, since 
the Gaussian wave function underestimates the peak den- 
sity iQ. In fact, the GP equation predicts that a soli- 
ton state for go = —2ir is stable (Fig. [l]). Although the 
Gaussian approximation does not accurately describe the 
exact soliton stability condition, Eqs. ( |l0| ) and ( pi] ) cap- 
ture the gi and J! dependences of i? m i n and u>hr- Fig- 
ure ^ (a) illustrates the monopole moment (r) versus 
(gi/Ci) 1 ' 2 , where the circles are obtained by varying g\ 
and the squares by varying ft. We note that the plots 
are fairly well proportional to (gi/fl) 1 ^ 2 , in agreement 



(a) 



t=0 it=0 

W 41 



(b) 



(c) 




t=0 




• 


• - 


-• 


• + 


t=13 
\ 




• 


• 



FIG. 3: Time evolution of the solitons with g(t) = — 2n + 
87rsin40i, where each soliton is prepared as in Fig. [j] (a). In 
(a) and (c), a weak radial harmonic confinement lo\ = O.lw 2 
is added. Two solitons are initially placed at x = ±3 with 
their phases differing by ir in (a) and with the same phase 
in (b). In (c), four solitons are initially placed on a slightly 
deformed square lattice, such that the adjacent solitons have 
phases differing by ir. The images are 8 x 4 in (a) and (b), 
and 8 x 8 in (c) in units of (ft/mwio) 1 ' 2 . 



with Eq. @. Figure || (b) shows the frequencies ui^ 
of the slow oscillations [as shown in Fig. [I] (b)] versus 
£l/gi- The plots are linear in £l/gi, which is consistent 
with Eq. (|ll|). 

The appropriate parameter range for stabilizing soli- 
tons with size R ~ 1 is found to be 0.4 < gi/il ^ 1.2 and 
go — — 2ir for f2 much larger than the breathing- mode 
frequency Wbr- For 17 < 10, the rapid oscillations are dis- 
turbed by nonadiabatic slow dynamics, and the soliton 
becomes unstable. In analogy with the inverted pendu- 
lum problem, this corresponds to a situation in which, 
the pendulum bob falls a long way during a single pivot 
cycle at a low pivot oscillation frequency. Thus, the effec- 
tive force fails to stabilize the pendulum bob. When <?i/Sl 
is small or \go\ is large, the effective force is not sufficient 
to prevent the atoms from accumulating at the center. 
As a result, the peak density first grows, the condensate 
then expands due to the R~ e potential, and subsequently 
most of the expanded atoms accumulate at the center due 
to the — i?~ 2 potential. Thus, the condensate repeatedly 
contracts and expands. In each expansion, atoms that 
are elastically scattered with high energy cannot return 
to the soliton region due to the absence of the external 
confinement potential, and the condensate gradually de- 
cays. 

We next consider the case in which two or more solitons 
coexist. The effective interaction between solitons de- 
pends on their relative phase . Figure || shows the dy- 
namics of solitons with g(t) = — 27r+87rsin40i. In Figs. H 



(a) and (c), a weak radial confinement w\ 
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added to keep the solitons from flying away. The initial 
states are taken to be ipsoiix + 3,y) =p ipaoi(x — 3,y) in 
Figs. |H (a) and (b), where ^> so ] is the wave function of a 
single soliton prepared as in Fig. [l| (a). When two soli- 
tons have phases differing by tt, they repel each other and 
oscillate at frequency ~ 2u>j_ [Fig. || (a)]. This situation 
is similar to the quasi- ID experiment performed by the 
Rice group ||, where adjacent solitons repel each other 
due to the 7r-phase difference. When two solitons have 
the same phase, they approach each other without the 
confinement potential, and merge into one condensate, 
which then expands [Fig. || (b)]. 

The repulsive interaction between solitons in 2D might 
suggest that a stable "soliton lattice" can be formed with- 
out an external periodic potential. However, this is not 
the case. Figure || (c) shows time evolution of a "soli- 
ton lattice," where initially the adjacent solitons have 
phases differing by tt and their configuration is slightly 
perturbed to a rhombus to examine the stability against 
lattice distortion. There is a dynamical instability such 
that diagonal solitons merge, which is followed by compli- 
cated dynamics. In the presence of a periodic potential, 
such alternative-phase structures are shown to be both 
stable in ID [Q and unstable, but can survive for a long 
time with appropriate parameters, in 2D and 3D ]l9| . In 
order to study the stability of our system in a periodic 
potential, we inserted the "plug" (an external potential 
such as e~ 4r ) between diagonal solitons, and found that 
the lattice structure survives for a long time. 

In conclusion, we have demonstrated that matter-wave 
bright solitons can be stabilized in 2D free space by oscil- 
lating the strength of the interaction around an attrac- 
tive value go < g cr with an amplitude g\ > \g \. The 
rapid oscillation of interaction produces an effective bar- 
rier that prevents the condensate from collapsing and sta- 
bilizes solitons. The interaction between solitons can be 
either repulsive or attractive, depending on their relative 
phase, and for repulsive interaction a soliton lattice can 
be stabilized using an optical plug. This novel technique 
of increasing dimensions of matter-wave bright solitons 
might be applied for quantum information processing us- 
ing BEC solitons on a microchip substrate pC[ |. It mer- 
its further study to examine whether BEC "droplets" in 
3D can be created with the appropriate parameters and 
ramp schemes. 
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